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Abstract 

This article presents new features of the MCSANC v.1.20 program, a 
Monte Carlo tool for calculation of the next-to-leading order electroweak and 
QCD corrections to various Standard Model processes. The extensions concern 
implementation of Drell-Yan-like processes and include a systematic treatment 
of the photon-induced contribution in proton-proton collisions and electroweak 
corrections beyond NLO approximation. There are also technical improvements 
such as calculation of the forward-backward asymmetry for the neutral current 
Drell-Yan process. The updated code is suitable for studies of the effects due to 
EW and QCD radiative corrections to Drell-Yan (and several other) processes 
at the LHC and for forthcoming high energy proton-proton colliders. 
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1 Introduction 

The forthcoming LHC data on Drell-Yan (DY) processes allows to access fi¬ 
nal states with very high invariant di-lepton masses, where the photon-induced 
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contributions become substantial relative to the standard quark-antiquark an¬ 
nihilation sub-processes. An accurate estimate of these contributions for hypo¬ 
thetical high mass resonance searches requires the inclusion into the theory pre¬ 
dictions sub-processes with photons in the initial pp state such as qj —)• 
qj qi~^£~ and 77 . 

Corrections to the neutral current Drell-Yan (NC DY) cross-section due to 
photon-induced process 77 —>■ i'^£~ can reach up to 10 — 20 % for high invari¬ 
ant mass Mg+i- with a choice of kinematic cuts typical for LHC experimental 
analysis. A first evidence of this kind of background was found by the ATLAS 
Collaboration in high mass NC DY cross section measurements [T]. 

Photon-induced Drell-Yan processes were carefully investigated in many 
papers, see for example [2] and [3]. 

In the TeV region of invariant masses the higher order two-loop electroweak 
(EW) and QCD leading terms are important for estimation of the theoretical 
uncertainties in different EW schemes. 

With this paper we continue the series of works dedicated to the devel¬ 

opment of MCSANC, a Monte Carlo tool based on the SANC modules [bj^ . 
We present an update of the integrator up to v. 1.20 with inclusion of the 
aforesaid corrections relevant for DY processes at the LHC at ^/s = 13 TeV. We 
briefly review the implementation into the framework of the MCSANC v. 1.20 
tool the following three new options: 

• photon-induced contributions. The implemented processes are: 

- q-f^ (for CC DY), 

- g 7 q£+£- (for NC DY), 

- 77 ^ £+£- (for NC DY); 

• Leading in two-loop EW and mixed EW( 8 )QCD radiative correc¬ 

tions; 

• forward-backward asymmetry 

This paper is organized as follows. In Section [2] we describe the implementa¬ 
tion of the photon-induced processes. Section [3] is devoted to the accounting of 
higher order radiative corrections using the p parameter at two loops. Results 
of our estimates of photon-induced processes, higher order radiative corrections, 
and forward-backward asymmetry A^Jp, as well as comparison between MC¬ 
SANC v. 1.20 and [ 2 ] and [7] are presented in Section [H A brief conclusion 
is given in Sect. [5l 

2 Photon-induced processes 

The introduction of photon-induced corrections into the SANC environment is 
presented in paper [3] where the MS subtraction scheme is realized. 

The DIS subtraction scheme can be realized using the sum of the following 

subtraction terms: _ _ 

(^MS^^piS^^(^MS^^DIS)^ ( 1 ) 

^This paper uses the same nomenclature as introduced in [3] 
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Here = 5i(ci), 5^® = (52 (ci) and the corresponding structure function 

are given by Eqs. ( 6 ) and ( 8 ) in paper [3]. The structure function 

Q^(l - 2 ;)/ 2 


= - log- 


(JiT 


xmf, 


-{x^ + (1 — xf) — 1 + 8x — 


( 2 ) 


should be used within the kinematics (5) of the above paper. Note that expres¬ 
sion ([1]) is for NC DY processes while only the first two terms in parentheses 
are present for CC DY ones. 

The 5^^® contribution is calculated within the following kinematics 

rv^ r-l-00 f-l r-l 

( 5 ^^^ = / dM / dY dxs / dcos-d qi [xi, fip)'y [x 2 , fip) 

g, J2mf J —00 Jo J—1 


d(T^T'(M,cos7?) 2M 


a 




dcos-d X 3 SQ 

v/so r-l-00 

dM / dY 
2mf J —00 


-yDIS 


0 (1 - xi) —- (5(1 - X 3 ) / C^,fiz)dz 


DIS/ 

-fli 


r-/: 


dcos-d 


2M 

•So 


qi {xi^dp) Q 0 (1 - xi)^ 


, 2 \ ^'^^^(-^jcos-d) Djg 

X 7 {x2,^iV) -773--(ics), 


dcos -d 

where the coefficient function reads 

T + x^ 


= - 


1 — X 


log 


1 — X 


+ 


9 + 5x 


and 


M = y^XiX2X3So, 


Y = - log-, 

2 X2 


M e 


Y 


M 


Xi = 


Xi = X1X3, X2 = 


-Y 


\/^ X 3 

where is the c.m.s. energy of colliding protons. Regularization of function 
Ca: at X —)• 1 is performed by plus-prescription explicitly shown in ([3]) . 

Another sub-process with photons in the initial state is 77 —)• i'^i~ which 
contributes to NC DY only. Here we give only the cross section of the process 
in the massless case, 


a 


'y'y^Born _ 


a^TT ft u 
u t 


and in the massive case, 


a 


77, 


,Born _ ^ 


1 




1-/3/ cos -d 
1 


1 + fii cos -d + 


4mj 


1 - 


1 -t- /3) 


1 — (di cos -d 


I + fdi cos -d 


1 — jdi cos -d + 


Arrij 


1 - 


l + (3f 


I + fdi cos -d 




( 4 ) 


( 5 ) 


mr 


where fdi = ^{s,mf ,mf) = y 1 — 4— and -d is the angle between the photon 
and outgoing lepton momenta in the center-of-mass system. 


3 


























3 Leading two-loop electroweak corrections 

In MCSANC v. 1.20 we follow the recipe introduced in Refs. [ 8 ] and [9], and 
later well described in Ref m- 

The p parameter is defined as the ratio of the neutral current to charged 
current amplitudes at zero momentum transfer, see for example [ 8 ]: 

^ Gnc{0) ^ 1 

Gcc{0) 1 — A/3 ’ 

where Gcc(O) = G^ is the Fermi constant dehned from the /r-decay width, and 
the quantity Ap is treated perturbatively 

Ap = Ap*-^^ + Ap^^^ + ... (7) 

Expanding ([ 6 ]) up to quadratic terms Ap^, we have 


p = 1 + Ap + Ap^ . 

The contribution to Ap, leading in G^m^ NLO EW, is explicitly given by 

Gu 


Ap 


( 1 ) 


= 


9 

mt 


( 8 ) 


(9) 


At the two-loop level, the quantity Ap contains two contributions: 

2as{ml) 


Ap = NcXt 1 -f p^^^ xt 


1 - 


97r 


■K + 3) 


( 10 ) 


They consist of the following: 

i) two-loop EW part at O(G^), second term in the hrst square brack¬ 
ets im, m and [9] with p^^^ given in Eq. (12) of [9] (actually, after discovery 
of the Higgs boson and determination of its mass, it has become sufficient to 
use the low Higgs mass asymptotic, Eq. (15), of [9]); 

ii) mixed two-loop EW( 8 )QCD at 0{G^as), the second term in the second 
square brackets, see in [H]- m for further details. 

From Eq. (j6]) , using intermediate vector boson propagators 1/(Q2 + M2), 
we derive: 


mz 


P = 


—9 9 ’ 

c^mj. 


( 11 ) 


where we have introduced a new parameter to distinguish from the usual 
for which we maintain the meaning /m^ to be valid to all perturbative 

orders. At the lowest order (LO) 

= 1 . ( 12 ) 




m 


w 




From Ea. (|lll) we have: 


c 2 = 


mz 


pm^. 


= (1 - Ap) cl 


(13) 


The universal higher order (h.o.) corrections, leading in may be 

taken into account via the following replacements: 




ac 


M „2 


s 2 - 


= si, + Ap c 


,2 

W 5 


= 1 - = (1 - Ap) C- 
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(14) 

(15) 














in the LO expression for NC DY cross section (see discussion after Eq.(3.48)) 
of [7]). 

As was argued in Refs. m,m and [ 8 ], this approach correctly reproduces 
terms up to O(Ap^). 

Given the replacements in Eq. (jl5[) . we get the following contributions of 
h.o. EW corrections to the scalar form factors d of the invariant amplitude. In 
the G^-scheme, 




vr 


(16) 


the form factor of 7 exchange effectively contains the factor 




acAtSw/'Sw = “G/i V (17) 

V / 


while four form factors of Z exchange, , contain a com¬ 

mon factor 1/c^. We consider as an example 




LL 


= «G, 


A„2 


T. 


■LL 


(18) 


Since the coupling factor aG^^js^ does not receive universal corrections, as 
follows from Ea. (ll7|) . we should insert only factors which come from 1/c^, 




1 -pLL 

■ A „2 „2 




^4^2 . 






4:S‘^ C 


(1 + Ap + Ap2). (19) 


In addition to Ea. (|19l) . form factors gf Z exchange contain the factor 

and form factor contains the factor s^. Therefore, form factors at 

NNLO order read: 

jzQQ ^ 1 + ^Ap, 

Sw 

= 1-FAp + Ap^ 

= (l + Ap + Ap2)(i + £|lAp), 

si 

j-QL _ (^1 _|__j__| ^Ap), 

si 

F^^^ = (1 + Ap +Api(l + %^Ap)i (20) 

st 


To avoid double counting one should remove the leading NLO EW contri¬ 
bution Q from the terms linear in Ap: Ap —;■ ^Ap — Ap^^^ ^ . 

Eas. (|20p were realized in new modules which compute the differential cross 
section contributions taking into account terms up to the order 0(Ap) (iho=l) 
or 0 (Ap 2 ) (iho= 2 ). 

We have shown analytically that the results obtained in this way agree with 
the corresponding expressions derived in [7]. 

^For the definition of the scalar form factors see Eqs.(29)-(30) in |16j . 
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4 Numerical results 


4.1 Phot on-induced processes 

In Tables [T] and [ 2 ] we present the LO inclusive cross sections for the process 

I ^ o t d'V 

pp ^ and the photon-induced contributions bgjn'y = - - for different 

(Jo 

ranges of the lepton pair transverse mass Mt,Iu and of the transverse momentum 
PT,i, respectively. For the sake of comparison we used the setup and input 
parameters from the paper [2]. Results of MCSANC (first line) are compared 
to the ones presented in Table 1 of [2] (second line). The lowest order cross 
sections are given in picobarns, the correction factors are shown in %. The 
numbers illustrate good agreement within the statistical errors of Monte Carlo 
integration. 


MT,iJGeV 

(To/pb 

(5q/q7/% 

50-OO 

4495.8(1) 

4495.7(2) 

0.047(3) 

0.052(1) 

100-cx) 

27.590(1) 

27.589(2) 

0.11(1) 

0.12(1) 

200-cx) 

1.7907(1) 

1.7906(1) 

0.24(1) 

0.25(1) 

500-CX) 

0.084696(1) 

0.084697(4) 

0.36(1) 

0.37(1) 

1000-cx) 

0.0065221(1) 

0.0065222(4) 

0.38(1) 

0.39(1) 

2000-cx) 

0.00027322(1) 

0.00027322(1) 

0.35(1) 

0.36(1) 


Table 1: Comparison of the LO and Sg/g^ for pp -P- in Mt^iu bins between 

MCSANC and [2]. 

In the same setup, in Tables [3] and 0] we present the LO inclusive cross 
sections for processes pp —)• e~i>eX and the photon-induced contributions 

In Table [5] we present the LO cross section for the process pp ^ e^e~X, 

(To in pb, and the corresponding contributions of photon-induced process 

(column 3) and 5^^ = (column 4). Here we used the setup and input 

o'o 

parameters given in [7]. The results of the MCSANC integrator are in the first 
rows, and the ones from Ref. [7] are in the second rows. Excellent agreement 
between these two calculations is observed. 

In the programs’ user interface these corrections are controlled by a new iph 
ffag added to the if lew parameter list. Setting iph =0 disables photon-induced 
contributions, iph=l includes the qy — ql'^l (NC DY) and (77 ^ (CC 

DY) components, and iph =2 includes the 77 ^ (NC DY) contributions. 

4.2 Higher order corrections 

In Table [5] we present the inclusive LO cross section gq in pb for the pro¬ 
cess pp e^e~X and the corresponding higher order corrections dh.o.weak = 


6 













PT,i/GeV 

(To/pb 

Sqlq^/% 

25-00 

4495.8(1) 

4495.7(2) 

0.059(3) 

0.065(1) 

50-OO 

27.590(1) 

27.589(2) 

4.6(1) 

4.7(1) 

100-cx) 

1.7907(1) 

1.7906(1) 

11.9(1) 

12.3(1) 

200 -cxo 

0.18129(1) 

0.18128(1) 

16.6(1) 

17.1(1) 

500-CX) 

0.0065221(1) 

0.0065222(4) 

16.2(1) 

16.7(1) 

1000-cx) 

0.00027322(1) 

0.00027322(1) 

13.1(1) 

13.5(1) 


Table 2: Comparison of the LO cross section and (5g/q7 for pp in px^i bins 

between MCSANC and [2]. 


MT,i-u/GeV 

(To/pb 


50-cx) 

3436.2(1) 

0.068(2) 

100-CX) 

20.037(1) 

0.113(4) 

200-CX) 

1.08169(1) 

0.223(2) 

500-CX) 

0.042127(1) 

0.328(3) 

1000-CX) 

0.002584(1) 

0.349(3) 

2000-cx) 

0.00008049(1) 

0.344(3) 


Table 3: The LO results and for pp ^ e in Mt^iu bins. 


Pt,i- / GeV 

(To/pb 

6,/,-,/% 

25-cx) 

3436.2(1) 

0.060(2) 

50-CX) 

20.037(1) 

5.30(1) 

100-CX) 

1.0812(1) 

16.22(2) 

200-CX) 

0.09503(1) 

26.31(2) 

500- CX) 

0.002584(1) 

34.87(4) 

1000-CX) 

0.00008049(1) 

39.72(2) 


Table 4: The LO results and 6q/q.^ for pp ^ e u^X in pt^i bins. 


o'h.o.weak 'pj^g setup and the input parameters are taken from 

o-Q 

paper [7]. The first rows represent results of the MCSANC integrator and 
the second ones show the numbers computed in [7]. Again we see excellent 
agreement between these two calculations. 

In the programs’ user interface the higher order corrections are controlled 
by a new iho flag added to the iflew parameter list. Setting iho =0 disables 


higher order correction contributions, iho=l includes linear 


(^Ap - Apb) 
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Mu/GeV 

(To/pb 

6qlq-y/% 

S-y-yfi/% 

Sh.o.weakf% 

50-OO 

738.813(5) 

738.773(6) 

-0.105(1) 

- 0.11 

0.17(1) 

0.17 

0.030(1) 

0.030 

100 -cxo 

32.7293(2) 

32.7268(3) 

-0.207(1) 

- 0.21 

1.16(1) 

1.15 

0.013(1) 

0.012 

200 -cx) 

1.48488(1) 

1.48492(1) 

0.381(1) 

0.38 

4.30(1) 

4.30 

-0.23(1) 

-0.23 

500-cx) 

0.080942(3) 

0.0809489(6) 

1.522(1) 

1.53 

4.92(1) 

4.92 

-0.29(1) 

-0.29 

1000-00 

0.0067998(1) 

0.00680008(3) 

1.901(1) 

1.91 

5.21(1) 

5.21 

-0.31(1) 

-0.31 

2000 -cxo 

0.00030375(1) 

0.000303767(1) 

2.343(1) 

2.34 

6.18(1) 

6.17 

-0.31(1) 

-0.32 


Table 5: The LO cross section (Tq oipp —)■ e^e~X in pb, the contributions of ( 57 ) and 
( 77 ) conhgurations in the initial pp state: 5q/q^ and and higher order corrections 
Sh.o.weak (iho= 2 ). Results of MCSANC (hrst rows) are compared to numbers from 
Ref. [7] (second rows). 


term contributions while iho=2 includes both linear and quadratic term 
contributions. 


4.3 Forward-backward asymmetry 

The forward-backward asymmetry is usually defined as (see my- 

F -B 


Afb = 


F + B 


where 


F = 


f 


da 


d cos d* 


-dcos d* 


B = 


i: 


da 


d cos d* 


-dcos d*. 


( 21 ) 


( 22 ) 


The cosine of the angle between the lepton and quark in the rest frame 
is then approximated by 

\Pz{l^l-)\ 2 


cos d* = 


X [p+(/-)p-(/+)-p-(r)p+(z+)] 


(23) 


where E is the energy, pz and pT are the longitudinal and transverse components 


1 


of the momentum vector, respectively, and p = —^{E Epz). 

V 2 

The MCSANC results for the forward-backward asymmetries (Figure [T]) 
were compared with the ones from Ref. m- A good agreement was found. 


5 Conclusion 

In this paper we have presented an update of the MCSANC integrator up to 
V. 1.20. The new features include systematic treatment of the photon-induced 





















Figure 1: The forward-backward asymmetry of pp —)■ X processes for the 

(left) and /i+/i“ (right) cases with bare event selection setups. 


contribution in proton-proton collisions and electroweak corrections beyond 
NLO approximation. The results of the calculations were compared with the 
results of other theoretical groups showing excellent agreement. The current 
version of MCSANC v. 1.20 is adjusted for studies of various effects due to 
EW and QCD radiative corrections to realistic LHC observables. 
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